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ABSTRACT 


This paper proposes the fractional-order proportional integral derivative 
(FOPID) controller, as a speed controller for permanent magnet direct current 
(PMDC) motor, instead of the traditional integer-order PID controller. The 
FOPID controller is the general form of the integer-order PID controllers, 
which found wide applications in all engineering fields. In this work a hybrid 
M-file and SIMULINK program is developed to simulate the overall system, 
the FOPID controller has five associated parameters. The optimum values of 
those parameters are found out by using particle swarm optimization 
technique. Simulation results show excellent command speeds tracking and 
superior dynamic response in conjunction with that of the integer-order PID 
controller. The proposed controller shows a high ability to overcome any 
external disturbance the system may be exposed; also, it performs a high 
degree of robustness to control the system in motoring and regenerative 
operating modes. 
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1. INTRODUCTION 

Fractional Order Proportional Integral Derivative controller (FOPID) is a most powerful type of PID 
controller in which the order of derivative and integral are fractional (PI*DY) rather than integer number. It 
offers an effective tool to descript the inheriting memory features of various articles. This is a great property 
of this controller in conjunction with ordinary PID controller; actually, this effect is not considered [1]. 

In addition, the fractional order controller is more convenient than integer-order models in the 
control of dynamical systems. Another advantage is the fact that fractional order controllers are insensitive to 
parameters variation of the system and also the controller parameters [2]. The FOPID control system presents 
a high degree of stiffness to overcome any disturbance or load variation, also its time and frequency 
responses can be adjusted efficiently as compared with the model predictive controller [3]. The fractional 
order controller (FOC) has already come in beneficial in many attractive engineering areas such as 
electronics, robotics, bioengineering, control theory and signal processing [4]. Concepcion A. M. et al, 
presented a self - tuning and auto-tuning to estimate the controller parameters of the FOPID controller, the 
method gave robust and more reliable control [5]. Rinu P. S. et al presented the Nelder-mead method and the 
ABC (artificial bee colony algorithm) to optimize the five parameters of Fractional-order PID controller [6]. 
Hassan N.A et al presented Pareto method for numerical optimization of the fractional order PID to give an 
estimation of the controller parameter to meet the engineering specification needs, and compared results by 
traditional PID controller [7]. 

From the literature, it can be noted that there are no previous attempts to investigate the non-linear 
behavior of the PMDC by parameters variation during operation. The resistance of the armature windings is 
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one of the main parameter whose has non-linear behavior due to the generated and rised heat in the windings 
by repetition of operation, which makes the system non-linear in nature. This paper proposes the use of the 
FOPID controller as a nonlinear controller of PMDC by taking into consideration the nonlinearty behavior of 
armature resistance, which is the main contribution of this work. This paper exploits the advanced features of 
the FO controller as an adjustable speed driver of a PMDC. The dominant parameters, of the proposed 
controller, are tuned by using the Particle Swarm Optimization (PSO) technique. The results of the proposed 
controller will be compared with that of the traditional integer-order PID controller. 


2. THE DIFFER-INTEGRAL AND FRACTIONAL ORDER SYSTEMS 

The fractional analysis is the sector that deals with the ability of powering the integral and 
differential operators by fraction number instead of an integer number. Mathematically, the fractional order 
integral and derivative can be defined in numerous ways. The general presentation of the differ-integrator 
operator 1s: 


ender) — “FO 
aD; f(t) = Gace (1) 





Where “gq” represents the real order of the differ-integral, “t” is variation domain, and “a” is the 
lower limit. In the Nineteen Sixties, Caputo presented a popular form to calculate the differ-integral. A 
fractional derivative of order “à” with variation in “t” domain and with initial condition t = O is become as 
[8]-[10]: 
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Where I'(z) is Euler’s Gamma function. When y < 0, then represent a fractional integral operator 
with —y order: 
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One wonderful benefit of using this definition is presenting easy explanations to initial and 
boundary conditions, in which the derivation of any constant value equals to zero. In the time domain, an n- 
term of fractional order differential equations are present the fractional order system: 


A, Dn y(t) +ay,_, D?r-1 y(t) + a, D” y(t) + ao D” y(t) = u(t) (4) 


Where D* = 2D? is the Caputo’s fractional derivation of order “A”. 
Tacking the Laplace transform to (4) the fractional order transfer function can be obtained, 
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Gn (s) = (5) 
Where; B, = (k=0,1,....,n)is a real number, 8, >8,-1>°:> 8, >8) >0 and a= 
(k = 0,1,....,n) is constant. 
The fractional derivative in Laplace form is: 


fr e-S* DY y(t)dt = SYY(s) — DR, SY-*1 y*(@) (6) 


The summation in the right-hand of (6) must be canceled if y < 0. 


3. FRACTIONAL ORDER CONTROL 

The PID controller is a universal controller based on a closed-loop mechanism has widespread 
applications in industrials controlling systems. The PID controller attempts to correct the error between an 
actual system feedback and a reference set-point by evaluating and then outputting a corrective action that 
could adjust the system according to that. The transfer function of an ordinary PID controller is: 
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The algorithm of the PID controller includes three independent parameters; the Proportion gain (K 
p) and the Integral and Derivative time constants (Ti & Td) respectively. K p evaluates the response to the 
instant error, while ( = ) estimates the integral reaction according to the summation of recent errors and 
1 
(Kp Ta) drive the derivative response according to the rate of changing errors. The algebraic summation of 
these actions is utilized to adjust the system by an appropriate controlling signal. 

In fact, the real controlled systems have fractional nature in general. Nevertheless, the fractionality 
could be small. Broadly, the approximate integer-order of the fractional system leads to a serious difference 
between the actual system and the mathematical model. Because of the fractional-order differential equations 
is a difficult calculus procedure then the integer-order models are used widely. The generalized form of the 
PID-controller, which also known as PI‘DY controllers, is as follow: 

Where, à and y are the integrator and differentiator orders respectively. 
Then, 


-r 
Gc (s) = K, (1+ — + Tas”), A, y > 0) (8) 
The time domain output is: 
u(t) = K, e(t) + K,D~e(t) + Ka DYe(t) (9) 


Where, K; = 2 and Ka = Kp Ta- 


In other words, the fractional order PI*DY controller is the global form of the conventional PID 
controller, which is a particular case controller. As shown in Figure 1, the PI*DY controller expands the 
conventional PID controller and extends it from point to surface. This extension increases the flexibility to 
controller adjustment, and performs more accurate response of the actual system. [8], [9]. 
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Figurel. Extend the PID control from point to surface 


4. MATHEMATICAL MODEL OF PMDC MOTOR 

Permanent magnet motors are used in various industrial applications because of their reliability, 
controllability, don’t need external excitation, and inexpensive [11], [12]. Permanent magnet DC motor is 
found in several low and high-power industrial applications. The equivalent circuit of the PMDC motor is 
shown in Figure 2, from which (10), (11) and (12) can be written. 
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Figure 2. The equivalent circuit of a permanent magnet DC machine 


Vin = laRa + La + KOO (10) 
T= Kyl; (11) 
JEZ = Te +T, + Bm Om (12) 


Where: Vin the supply voltage (V), Ia current of the armature (A), @ is the pole magnetic flux (Wb), 
(Ra) resistance of the armature (Q), (La) inductance of the armature (H), K machine constant (V.s/rad), om 
rotor speed (rad/s), Te developed torque(Nm), TL Load torque (Nm), Bm damping coefficient (N.m.s/rad), 
and J inertia constant (kg.m2). Clearly, from (10), the armature resistance and inductance are the dominant 
parameters which judged the transient and steady state response. The variation of armature resistance is taken 
periodically about 10% of its cold value as shown in Figure 3. 

















Figure 3. The variation of armature resistance 


5. PARTICLE SWARM OPTIMIZATION (PSO) 

The particle swarm optimization is firstly presented by James and Russell in 1995[13], a seeking 
algorithm having a random probability distribution or pattern inspiring by a community habit of the flock of 
birds [14], it is one of the powerful optimization techniques. It has been observed to be strong in nonlinearly 
optimization issues [15], [16]. An improved PSO was presented in 1998 to upgrade the performance of the 
conventional PSO. An addition coefficient known as inertia weight is introduced [14]. PSO is a set of rules 
associates both cognition and public collaboration of birds. 

In this method, each bird "particle" searches for the optimum position by updates its location within 
the swarm, according to the advantage of its memory of the best position and knowing about the global best 
location [15]. The optimization procedure obtained by (13) and (14). 


Vi (k + 1) =W vi(k) T Cı Ry (Ibest — xi(k)) + C2 R2 (Poest — Xi (k)) (13) 
xi(k +1) = x,(k) + v;(k + 1) (14) 
i=1,2,...,n 


Where: vi is the ith particle velocity, xi is the ith particle position, k iteration number, C1, C2 known 
as the cognitive and social coefficients, w inertia weight factor, Rl, R2 random variables of from 0 to 1, 
Popesti dividual best position of particle 1, Zpest best global position of all the particles in the swarm, and n 
number of birds (particles). 

If the condition in (15) is met, then the position is updated through (16): 


f Xix ) < f (Prest) (15) 
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Where: f perform the minimization objective fitness function. The flowchart of the PSO algorithm is 
illustrated in Figure 4. 

Usually, in optimization techniques, there are several fitness criteria to evaluate system performance 
as such as Integral Absolute Errors (IAE), the Integral of Square Errors (ISE), and Integral of Time Square 
Errors (ITSE). These performance criteria are inclusive the overshoot, rising time, settled time and steady- 
state error. In addition, it has been specified the optimization process and the tightness of the driving system 
[17]. In this work the ITSE fitness function, illustrated in (17), is used as the performance criteria of the 
output response of the system: 


ITSE = f. te?(t).dt (17) 


The block diagram of the optimization process for FOPID-PSO system is shown in Figure 5. A 
MATLAB program is developed to find out the optimal values of the five controller parameters à, y, Kp, Ti, 
and Tg. A minimization algorithm is applied in the search domain of the particles position and velocity, the 
optimum values of the FOPID controller are obtained in 100 iterations as shown in Figure 6. 


Initial random positions and velocities of the 
particles 
Evaluate the fitness of each particle 
Obtain best particle’s position 


Update the global best position over all 
previous best positions 


Evaluate the new position 


ls the no. of iterations reaches 





Figure 4. Flow chart of the PSO algorithm 
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Figure 5. The block diagram of the optimized FOPID controller 
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Figure 6. The controller parameters during the optimization process 


6. SIMULATION RESULTS 

The overall system simulation is shown in Figure 7. A MATLAB-SIMULINK model is designed; a 
built-in M-file is developed to calculate the fractional order integral and derivative in a specific sampling 
time. The speed, torque, and armature current response to rated speed command of the PMDC motor under 
the full-load condition are depicted in Figure 8. Also, the response of the system to a stair steps input of 
different speed commands under full-load condition can be shown in Figure 9. 
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Figure 7. The overall system simulation 
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Figure 9. Stair response of the FOPID controller 


In order to realize the features of the proposed controller, a comparison performance is investigated 
in conjunction with the traditional integer-order PID controller. The performance for both rated of speed and 
stair commands are shown in Figure 10 and Figure 11. Moreover, the robustness of the proposed controller 
can be proven in four-quadrant operation conditions of the PMDC motor. Positive and negative speed 
commands are applied to the motor under various positive and negative load torques. The speed, torque, and 
current performance of the system for this operation are shown in Figure 12. 
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Figurel 1. A comparison performance of different speed command 
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Figure 12. Four-quadrant operation 
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7. DISCUSSION 

Obviously, the performance of the PMDC under the action of the FOPID controller shows an 
excellent response, for both of rated step and stair steps input commands, although the non-linear behavior of 
the msystem. Figure 8 and Figure 9 show accurate speed tracking for at no-load and full-load conditions, the 
performance shows a superior dynamic response. Also, the comparison investigation, between the fractional- 
order and the integer-order PID controllers, shows eminent features of the proposed controller. Figure 10 and 
Figure 11 illustrate the dynamic response (maximum overshoot, rise time, settling time, and steady-state 
error) to the rated step input of the integer-order PID controller about 26%, 0.23s, 1s, and 1.3% respectively. 
Whereas, the dynamic responses of the FOPID controller are: 1.7%, 0.23s, 0.37s, and 0.2% respectively. In 
addition, the starting current period of the FOPID is about half of that of the integer-order PID, this will 
reduce the winding temperature due to multi-starting operation. Moreover, the transient response of the 
traditional PID controller due to full-load disturbance shows: maximum overshoot = 2.7% and settling time = 
0.5s. Whilst, the FOPID controller shows: maximum overshoot = 0.7% and settling time = 0.25s. Finally, the 
four-quadrant operation of the PMDC motor in Figure 12 shows a high degree of robustness for different 
positive and negative speed and load torque. A superior speed tracking is performed under both of delivered 
and consumed power modes. 


8. CONCLUSION 

This paper proposed the FOPID controller, instead of the integer-order PID controller, to control the 
speed of the PMDC with considering the non-linear behavior of armature resistance. Both of the fractional- 
order and the integer-order PID controllers are simulated, and the responses of them are investigated under 
different operating conditions. Simulation results show excellent command speeds tracking and superior 
dynamic response although the variation of system parameter. The proposed controller shows a high ability 
to overcome any internal or external disturbance the system may be exposed; also, it performs a high degree 
of robustness to control the system in motoring and regenerative operating modes. 
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